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1 Homogeneous media 

Suppose we have Poisson's equation with constant coefRcients (= 1 for simplicity), i.e. in a 
homogeneous medium, i.e. in "empty space", i.e. 

in some region fi, with boTindary conditions so that this has a unique solution, e.g. Dirichlct 
boundaries u\dn = 0. Then we can write the solution in terms of the Green's function 
Go(x,x'): 

«(x)= / Go(x,x')/(x')d"x', (1) 

in n dimensions, where — V^Go(x, x') = 5(x — x'). For example, if = M'^, then Go(x, x') = 
1 

47r|x— x'l ■ 

2 Inhomogeneous media 

Now, suppose we have non-constant coefRcients c(x) > 0, corresponding to an inhomoge- 
neous medium (i.e. different "materials" at different points in space). This could enter 
Poisson's equation in several ways, for example: 

1. — cV^w = / — for example, in a stretched string or drum, 1/c would be proportional 
to the density, so c(x) would represent a variable density. 

2. — V-(cVm) = / — for example, in electrostatics ^/c would be proportional to the refrac- 
tive index; in a stretched string or drum c would be proportional to a tension/elasticity; 
in a diffusion problem c would be a diffusion coefficient; in a heat-conduction problem 
c would be proportional to a thermal conductivity. 

3. —V^u+cu = f — for example, in quantum mechanics this would represent Schrodinger's 
equation with a variable potential energy c. 

4. Various other generalizations and combinations, e.g. a multidimensional Sturm- 
Liouville equation, — ciV • (c2Vu) + c^u = f, for different functions Ci_2,3(x) (and 
C2 could even be a positive-definite matrix). 

All of these are of the form Au = f where A is self-adjoint and positive-definite (assuming 
zero Dirichlet boundary conditions and an appropriate inner product {u,v)), as we've seen 
previously in class, so they should have unique solutions (excluding pathological c functions). 
However, the solutions may in general be quite different from those of — V^u = /. Can we 
relate them to Go, the Green's function for empty space? 
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2.1 -cV^u = / 

This case is quite trivial to relate to empty space: just rewrite it as — V^u = //c (valid since 
c > 0), in which case we can use eq. (1) to obtain 

«(x)= [ Go(x,x')^d"x', 

from which it follows that the Green's function of this problem is G(x,x') = Go(x,x')/c(x'). 

2.2 -V • (cVu) = f 

To make it look like empty space, we employ the product rule to write V(cVu) = (Vc) • 
(Vu) + cV'^u, obtaining: 

V72 / Vc 
—V u= — I • yu. 

c c 

If we substitute the rigiit-liaud side as if it were "/" in the empty-space problem, from 
eq. (1), and note that (Vc)/c = Vine, we obtain: 

m(x)= / Go(x,x') 
Jn 

Notice, however, that the unknown u appears on both the left- and right-hand sides. This is 
a volume integral equation^ for m(x). There are various numerical methods to solve such 
"VIE" problems approximately by discretizing space, but we won't cover these in 18.303. 
However, there are still several lessons to be learned from this equation. First, one can think 
of the inhomogeneous solution as the sum of "homogeneous" solutions Go(x,x') from the 
right-hand-sides /(x')/c at x' plus a "scattered" solution from the solution u creating 
new source terms at inhomogeneities Vc. Second, there are various general situations 
where this VIE simplifies considerably. In the following, denote 

tio(x) = [ Go(x,x')4^d"x', 

the part of m(x) that doesn't come from the inhomogeneity Vc, the "incident" solution. 
We then can write u = uq + Bu where B is the linear integral operator corresponding to 
Bu = ^ GoV'lnc- V'w. 

(This terminology of "incident" and "scattered" parts of the solution has its origins in 
wave-equation problems, where the solution represents a wave propagating outwards from a 
source /, and is actually bouncing off of objects/inhomogeneities. In the present problem, 
there is no time dependence so this terminology is only an analogy.) 

2.2.1 Piecewise-homogeneous media 

The most common and important case of inhomogeneous media is where c(x) is piecewise- 
constant. e.g. in one region you have glass, in another region you have metal, and in 
another region you have water, each corresponding to a different constant value of c. 

For example, suppose we have = M^, with —V - (cV?i) = / for the situation depicted in 
figure 1: c(x) = Ci in some volume V and = C2 outside of V. In this case Vine is actually 
a delta function at the interface, multiplied by the magnitude lnc2 — Inei = ln(c2/ci) of 

'^More specifically, it is a "second-kind" integral equation, since u appears both inside and outside the 
integral. 



c(x') 



-V'[lne(x')] • V'w(x') 



rf"x'. 



(2) 
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Figure 1: Schematic example of piecewise-constant coefficients c(x): c(x) = ci in some 
region V (with boundary dV, and n denoting the outward unit-normal vector) , and c(x) — C2 
otherwise. 



the discontinuity, also multiplied by a unit-normal vector n at each point on the interface 
(giving the direction Vine). Then the VIE (2) simplifies to uq plus a surface integral: 

u{k) = uo(x) + ln(c2/ci) S Go(x, x') V'u(x') • dA' , (3) 

JJdV 

where dA' = hdA' is the usual outward-normal differential area in a surface integral. This 
is now a surface integral equation (SIE) for m(x): once we know n- Vu on the surface, we 
can get u(x) everywhere! Physically, this can be interpreted as "scattering" of the solution 
off the interface between ci and C2, as represented by "source" terms at every x' £ dV. 

However, this isn't quite right. The problem is that Vw • h isn't actually continuous 
across the interface in general, and so evaluating it on the surface dV is not well-defined. 
We can see this by looking at the original equation V • (cVn) — f: unless / happens to have 
delta functions right on the dV interface, the quantity cVm must be continuous across the 
interface in the h direction. That is, cVu • h is continuous, not Vw • h. So, before we can 
evaluate a surface integral, we need to rewrite things in terms of cVu • h. We can do this, 
since 

VRncl • V-u = — • Vii ^ • cVu = V f -- ) • cVm 
c c-' V c / 

in eq. (2). And V(— 1/c) is also a delta function multiplied by n and the magnitude ^ ^ ^ 
of the discontinuity, and hence we obtain a corrected SIE: 

u{^) = ^.o(x) + - ^) Go(x, x') [c(x') V'u(x')] • dA', (4) 

which is well-defined since the integrand is now a continuous quantity across the interface. 
(And even this is not quite in the form that is used in numerics.) 

Like our SIE equations for the case where itself has a boundary, these SIE equations are 
solved by parameterizing the surface unknowns via some discretization, and then choosing 
the unknowns so that u satisfies appropriate continuity conditions at the interface dV. The 
details of this get rather tricky very quickly, but these can yield very efficient computational 
methods because they only involve unknowns at interfaces and handle homogeneous regions 
(even infinite homogeneous regions) analytically via Gq. 
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2.2.2 Born approximation 



We have written the solution in the form u = uq + Bu, where uq is the solution due to / 
ignoring the inhomogeneity, and B is an integral operator giving the "scattered" portion of 
the solution Bu due to the inhomogeneity. We can then formally write: 

u = uo + Bu = uo + B{uo + Bu) = uo + Buq + B^{uo + Buq + B^u) 

\k=0 / 

which is called a "Born-Dyson" series/expansion (sometimes omitting one name or the 
other). Equivalently, (1 — B)u = uq implies u = (1 — B)~^uo, and the Taylor series for 
(l-S)-i is EkB". 

What does this represent physically? uq is the "incident" portion of the solution, before 
scattering off the inhomogeneity. IBuq is the portion of the solution where this "incident" 
solution has scattered once off the inhomogenity, producing a scattered solution 

Buo = [ Go(x,x')V'lnc(x')- V'uo(x')rf"x'. 

But of course, this portion of the solution also scatters off of the inhomogeneity, producing 
a portion of the solution representing incident solutions that have scattered twice: 

B^uo= [ Go(x,x')V'lnc(x') • V 
Jn 

i.e. a 7/0 at x" produces a source term due to V'c, which "travels" from x" to x' via 
G'o(x', x"), then produces a source at x' due to V'c, then "travels" from x' to x via Go(x, x'), 
and of course this must be integrated over all possible scattering points x' and x". Then 
B^uq represents things scattering three times, and so on.^ 

Now, suppose we have a system that is nearly homogeneous, e.g. Vc is small. Then one 
might expect that the scattered portion Bu of the solution to be small. In this case, we may 
be able to approximate this series by keeping only the first two terms — if scattering once 
has small amplitude, then scattering twice should have even smaller amplitude, (e.g. if Vc 
is small, this could be thought of as expanding as a power series in Vc, since B'^ ^ \Vc\''.) 
This is the Born approximation: 

u(x) pa Uo(x) + Buo, 

and is an extremely useful way to think about nearly homogeneous problems. 



/ 



Go(x',x")V"lnc(x") • V"uo(x")d"x" 



3 Example: Inhomogeneity in a small volume 

Suppose we are solving —V • (cVm) = / in Q = with a point source /(x) = d{x — Xq) 
at xq. Furthermore, suppose that c(x) is piecewise-constant as in figure 1, with c(x) = C2 
everywhere except in a volume V, centered at xi, where c(x) = ci. Now, suppose that we 
want the solution m(x), but are far from V: both the source point xq and the desired point 
X are far from V, with |xi — xo| and |xi — x| both much bigger than the diameter of V. This 
is shown schematically in figure 2. In this case, we should expect the effect of the "scattered" 

^In quantum mechanics, this kind of series of events is sometimes represented graphically by the notation 
of "Feynman diagrams," and can be generalized to nonlinear problems and other effectivei inhomogeneities. 
In that context, this process of summing all possible scattering sequences is sometimes mysteriously described 
as a "particle exploring all possible paths between two points," but is really just a consequence of particles 
being described by PDEs (Schrodinger's equation, in single-particle quantum mechanics), rather than ODEs. 
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Figure 2: Schematic of problem with an inhomogeneity in a small volume V (centered at 
Xi): we have a source at Xq and want the solution at x, with both Xq and x much farther 
from Xi than the diameter of V. 



solution from V to be small at x, and a Born approximation should apply. Furthermore, 
we will assume ci ~ C2 (though not exactly equal!), so that we can neglect the effect of 
the discontinuity in Vu mentioned after equation (3) above (which greatly complicates the 
application of any Born-like approximation in this problem because it would prevent us from 
using u uq in V).^ 
In this case, 

uo(x) = Go(x,xo)/c(xo) = ■ r, 

47rc2|x - Xo| 

so in the Born approximation we write: 

u(x) « Uo(x) + Buq, 

where the scattered part of the solution, applying the SIE form (4) [valid when ci ~ C2], is 



Buq = ln(c2/ci) 
= ln(c2/ci) 
= ln(c2/ci) 



Go(x,x')V'uo(x') -dA' 



dV 



V • [Go(x,x')V'wo(x')]rf^x' 
V'Go(x,x') •V'uo(x') + Go:5 



where in the second line we applied the divergence theorem, and in the third line the product 
rule led to a V^Uq term, where V^Uq — —S{x — Xg) is zero in V (since Xq is outside of V). 

Now, since V is small compared to the distance from x and Xq, the distances |x' — x| and 
|x' — xol hardly change for any x' G V, and so the V'Gq and V'mq terms are approximately 
constant in this integral and we can just pull them out, giving the approximation: 

Buo « ln(c2/ci) V'Go(x,x') • V'uo(x')lx'=xi volumc(y). 
We can compute these gradients explicitly: 

1 x'-y 



|x' -y| 



|x' -y| 



and hence: 



u(x) 



47r X — Xi 



01 



+ ln(c2/ci) 



(Xl - X) ^ (Xl - Xq) 

47r|xi-xP 47rc2|xi-xo| 



■ volume (F). 



(5) 



^It turns out that many people get this wrong in electromagnetism for cases when ci and C2 are very 
different, as discussed in my paper on a closely related subject, "Roughness losses and volume-current 



methods in photonic-crystal waveguides,' 
"stevenj /paper s/JolmsonPo05 . pdf 



Appl. Phys. B 81, 238-293 (2005): http://math.mit.edu/ 
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Notice that the ampHtude of the scattered term vanishes as volume(y) 0, as expected. 
Notice that it also depends on the sign of (xi — x) • (xi — xo). Why is that? What does a 
V'Go source "mean," physically? 

3.1 Dipole sources 

Consider the following problem in = M'^, requiring as usual that solutions vanish at oo: 

~V^Dp{Ti, x') = -p • V5(x - x') = +p • V'(5(x - x'). 

This is like the Green's function equation, except now we have put the derivative of a delta 
function on the right-hand side, with some constant vector p (the "dipole moment"). Recall 
what the derivative of a delta function is: 

hp . V<5(x - ^')m = [5(x - x')]{p • V0} = p . = ^^^ H^' + ep)-^<l>{^'-ep) ^ 

and hence (similar to pset 5), 

, (?(x - x^ - ep) - (5(x - x' + ep) 
— p • Vdix — X = lim . 

That is, the derivative of a delta function is a limit of limit of two delta functions of opposite 
sign, displaced proportional to p. In 8.02, where delta functions are "point charges," this is 
what you would have called an "electric dipole." 

We can solve for Dp quite easily, because we know the solution Gq to — V'^Go(x,x') = 
(5(x — x'), and V and V' derivatives can be interchanged in their order: 

-p • V,5(x - x') = p • V [<5(x - x')] = p • V [-V2Go(x, x')] = -V^ [p • V'Go(x, x')] , 

and hence ^ 
i?p(x, x') - p • V'Go(x, x') = p • 

47r|x — x'l'' 

In electrostatics, ths would be the potential of a dipole. Note that this falls off as ^ 
l/|x — x'p, whereas Gq falls off as ~ l/|x — x'|. 

Given this solution, we can now interpret the scattered part of the solution (5) above: a 
small inhomogeneity gives an effective dipole source p at Xi, where 

p = -ln(c2/ci)— ^^-^ — volume(V^). 

In electrostatics, for a typical case where 1^ is a small piece of matter in vacuum, C2 < ci, 
so p is parallel to xi — xq. Physically, a positive point charge induces a dipole moment p 
pointed away from the charge, because a "+" charge at xq pushes "+" charges in V away 
from it, as shown below. 

+/ 
+/ 
+/ 
+/ 

V /- + 
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